FIX: x0

#ReactionName
#reaction
#reactionRate

R1:
s1 + s2 = s3
k1*s1*s2 - k2*s3

R2:
s3 = s4
k3*s3 - k4*s4

R3:
s4 = s1 + s5
k5*s4 - k6*s1*s5

R4:
s5 + s6 = s7
k7*s5*s6 - k8*s5-s6

R5:
s7 =s8
k9*s7 - k10*s8

R6:
s8 = s2 + s9
k11*s8 - k12*s2*s9

R7:
s9 + s10 = s11
k13*s9*s10 - k14*s11

R8:
s11 = s12
k15*s11 - k16*s12

R9:
s12 = s13 + s6
k17*s12 - k18*s6*s13

R10:
s13 + s14 = s15
k19*s13*s14 - k20*s15

R11:
s15 = s16
k21*s15* - k22*s16

R12:
s16 = s13 + s17
k23*s16 - k24*s13*s17

R13:
s17 + s13 = s18
k25*s17*s13 - k26*s18

R14:
s18 = s19
k27*s18 - k28*s19

R15:
s19 = s17 + s10
k29*s19 - k30*s17*s10

R16:
s13 + s20 = s21
k31*s13*s20 - k32*s21

R17:
s21 = s22
k33*s21* - k34*s22

R18:
s22 = s13 + s23
k35*s22 - k36*s13*s23

R19:
s23 + s9 = s24
k37*s23*s9 - k38*s24

R20:
s24 = s25
k39*s24 - k40*s25

R21:
s25 = s23 + s6
k41*s25 - k42*s23*s6

R22:
x0 = s1
k43*x0 - k44*s1


k1 = 10.0000
k3 = 0.0100
k5 = 1.0000
k7 = 10.0000
k9 = 0.0100
k11 = 1.0000
k13 = 8.4770
k15 = 14600.0000
k17 = 6.7000
k19 = 10.0000
k21 = 0.0010
k23 = 1.0000
k25 = 10.0000
k27 = 0.0140
k29 = 1.0000
k31 = 10.0000
k33 = 0.0010
k35 = 1.0000
k37 = 10.0000
k39 = 0.0100
k41 = 1.0000
k43 = 1.0e-16

k2 = 1.0000
k4 = 0.0100
k6 = 10.0000
k8 = 1.0000
k10 = 0.0100
k12 = 10.0000
k14 = 0.3325
k16 = 16.0000
k18 = 0.0600
k20 = 1.0000
k22 = 0.0001
k24 = 10.0000
k26 = 1.0000
k28 = 0.0001
k30 = 10.0000
k32 = 1.0000
k34 = 0.0001
k36 = 10.0000
k38 = 1.0000
k40 = 0.0001
k42 = 10.0000
k44 = 1.0e-16

x0 = 1.0

s1 = 0.01 
s2 = 0.1
s3 = 1e-4
s4 = 1e-4
s5 = 1e-4
s6 = 1.4
s7 = 1e-4
s8 = 1e-4
s9 = 1e-4
s10 = 1.0
s11 = 1e-4
s12 = 1e-4
s13 = 1e-4
s14 =0.1
s15 = 1e-4
s16 = 1e-4
s17 = 1e-4
s18 = 1e-4
s19 = 1e-4
s20 = 0.1
s21 = 1e-4
s22 = 1e-4
s23 = 1e-4
s24 = 1e-4
s25 = 1e-4
